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ABSTRACT 

We use the angular power spectrum of the Cosmic Microwave Background, measured during the North 
American test flight of the Boomerang experiment, to constrain the geometry of the universe. Within 
the class of Cold Dark Matter models, we find that the overall fractional energy density of the universe, 
n, is constrained to be 0.85 < 12 < 1.25 at the 68% confidence level. Combined with the COBE 
measurement and the high redshift supernovae data we obtain new constraints on the fractional matter 
density and the cosmological constant. 

Subject headings: cosmology: Cosmic Microwave Background, anisotropy, measurements, power 
spectrum 



1. INTRODUCTION 

The dramatic improvement in the quality of astronom- 
ical data in the past few years has presented cosmologists 
with the possibility of measuring the large scale prop- 
erties of our universe with unpre cedented precision (e.g. 
Kamionkowski & Kosowsky 1999). The sensitivity of the 
angular power spectrum of the Cosmic Microwave Back- 
ground (CMB) to cosmological parameters has lead to 
analyses of existing datasets with increasing sophistica- 
tion in an attempt to measure such fundamental quantities 
as the energy density of the universe and the cosmologi- 
cal constant. This activity has lead to improved methods 
of Maximum Li kelihood Estimation (Bond, Jaffe & Knox 



most recent data: again, the 68% likelihood contours com- 
fortably encompass the Einsten-de Sitter Universe. All 
these previous analyses were restricted to the class of open 
and flat models. 

In this letter we present further evidence for a flat uni- 
verse from the CMB. Usi ng the methods for parame- 
ter estimation described in BJK98, we perform a search 



1998, hereafter pjK98 
enlar g ing the range o 



Bartlctt et al 1999| ), attempts at 
possible parameters ( Lineweaver 



1998| , [Tegmark 1999| , [Melchiorri et al. 1999| ), and the 
incorporation of systematic uncertainti es in the experi- 
ments (D odelson & Knox, from now on DK99, Ganga et 



al . 1997D . 

Within the class of adiabatic inflationary models there 
is now strong evidence from the CMB that the universe is 
flat. The most extensive range of parameters has been con- 
sidered by Tegmark (1999) where the author found that 
a flat universe was consistent with CMB data at the 68% 
co nfidenc e level. A more thorough analysis was performed 
incorporating the non-Gaussianity of the likeli- 



in cosmological parameter space for the allowed range 
of values for the fractional density of matter, Hm, and 
cosmological constant, J7a, given the recent estimate of 
the angular power spectrum from the 1997 test flight of 
the Boomerang ex periment (see the companion paper 
Mauskopf et al. 1999| ). We obtain our primary constraints 
from this data set alone and find compelling evidence, 
within the family of adiabatic inflationary models for a flat 
universe. In section]^ wc briefly describe the Boomerang 
experiment, the data analysis undertaken and the charac- 
teristics of the angular power spectrum obtained. In sec- 
tion H we spell out the parameter space we have explored, 
the method we use and the constraints we obtain on the 
fractional energy density of the universe, = ^Im + ^^A- 
Finally in section | we discuss our findings and combine 
it with other cosmological data to obtain a new constraint 
on the cosmological constant. 
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hood function, possible calibration uncertainties and the 
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2. THE DATA 

The data we use here are from a North American 
test flight of Boomerang (Boomerang/NA), a baUoon- 
borne telescope designed to map CMB anisotropics from 
a long-duration, balloon-borne (LDB) flight above the 
Antarctic. A detailed descrip tion of the instrument can 
be found in Masi et al. 1999. A description of the data 
and observations, with a discussion of calibrations, sys- 
temat flc effects and sig nal reconstruction can be found in 



Maus kopf ot al. 1990 . This tost flight produced maps of 



the CMB with more than 200 square degrees of sky cover- 
age at frequencies of 90 and 150 GHz with resolutions of 
26 arcmins FWHM and 16.6 arcmins FWHM respectively. 

The size of the Boomerang/NA 150 GHz map (23,561, 
6' pixels) required new methods of analysis able to incor- 
porate the effects of correlated noise and new implementa- 



tions 



capable of processing large data sets. The pixclizcd 



map a nd angular power spectrum were produced using the 
MADCAP software package of BorriU (1999a, 1999b) (sec 



http://cfpa.berkeley.edu/'-~^borrill/cmb/madcap.html) on 
the Cray T3E-900 at NERSC and the Cray T3E-1200 at 
CINECA. 
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Fig. 1. — The Likelihood function of each of the eight band 
powers, 4fl = (58, 102, 153, 204, 255, 305, 403, 729), reported in 



ansatz of BJK9J 



Mauskopf et al. 199S: computed using the offset lognormal 



The angular power spectrum, C^, resulting from the 
analysis of the 150 GHz map was estimated in eight bins 
spanning t with seven bin's centered between ^ = 50 to 
i = 400 and one bin at ^ = 800. 

The bin correlation matrix is diagonalized as in BJK98 
resulting in eight orthogonalizcd (independent) bins. We 
present the likelihood for each orthogonalizcd band power 
in Figu re 1, using the offs et lognormal ansatz pr oposed in 
BJK98i As described in [Mauskopf et al. 1999| the data 



show strong evidence for an acoustic peak with an ampli- 



tude of ~ 70/zKcAfB centered a.t I ^ 200. 

3. MEASURING CURVATURE 

The Boomerang/NA angular power spectrum covers 
a range of £ corresponding to the horizon size at decou- 



pling. The amplitude and shape of the power spectrum 
is primar ily sensitive to the overall curvature of t he uni- 
verse, r2 ( Doroskevich, Zeldovich, & Sunyaev 1978); other 
parameters such as the scalar spectral index, ns, the frac- 
tional energy density in baryons, O^, the cosmological 
constant, ^a, and the Hubble constant, Hq = lOOh km 
sec~^, will also affect the height of the peak and there- 
fore some "cosmic confusion" will arise if we att empt in- 
dividua l constraints on each of the parameters (Bond et 
al 1994). In our analysis we shall restrict ourselves to the 
family of adiabatic, CDM models. This involves consid- 
erable theoretical predjudice in the set of parameters we 
choose to vary although, as the presence of an acoustic 
peak at £ ~ 200 becomes more certain, the assumption 
that structure was seeded by primo rdial adiabat ic pertur- 
bations becomes more c ompelling ([Liddle 1995; however 



counterexamples ex ist, Turok 1997| , purrer fc Sakellari- 
adou 1997| , |Hu 1999| ) 

We should, in principle, consider an 11-dimensional 
space of parameters; sensible priors due to previous con- 
straints and the spectral coverage of the Boomerang/NA 
angular power spectrum reduce the space to 6 dimensions. 
In particular, we assume Tc — (lacking convincing evi- 
dence for high redshift reionization), we assume a neglige- 
able contribution of gravitational waves (as predicted in 
the standard scenario), and we discard the weak effect due 
to massive neutrinos. The remaining parameters to vary 
are ftcDM, ^A, ^b, h, ns and the amplitude of ffuctu- 
ations, Cio, in units of CyP^^ . The combination fls/i^ 
is constrained by primordial nucleosynthesis arguments: 
0.013 < Q.Bh'^ < 0.025, while we set 0.5 < h < 0.8. For 
the spectral index of the primordial scalar ffuctuations we 
make the choice 0.8 < ns < 1.3 and we let a 20% variation 
in CiQ. As our main goal is to obtain constraints in the 
— QcDM + ^B, ^a) plane, we let these parameters 
vary in the range [0.05, 2] x [0, 1 



Proceeding as in DK99 



we attribute a likelihood to a point on this plane by find- 
ing the remaining four, "nuisance" , parameters that maxi- 
mize it. The reasons for applying this method are twofold. 
First, if the likelihood were a multivariate Gaussian in all 
the parameters, maximizing with regards to the nuisance 
parameters corresponds to marginalizing over them. Sec- 
ond, if we define our 68%, 95% and 99% contours where 
the likelihood falls to 0.32, 0.05 and 0.01 of its peak value 
(as would be the case for a two dimensional Multivariate 
Gaussian), then the constraints we obtain are conserva- 
tive relative to any other hypersurface we may choose in 
parameter space in the sense that they rule out a smaller 
range of parameter space than other usual choices. 

The likelihood function for the estimated band powers 
is non-Gaussian but one c an apply the "radical compres- 



sion" method proposed by BJK98; the likelihood function 



is well approximated by an offset lognormal distribution 
whose parameters can be easily calculated from the output 
of MA DCAP. The theory CjS are generated using CMB- 
FAST ( Beljak fc Zaldarriaga 1996| ) an d the recent imple- 
mentation for closed models CAMB ([Lewis, Challinor fc 



Lasenby 199£). We search for the maximum along a 4 di- 



mensional grid of models, using the fact that variations in 
Cio and Us are less CPU time consuming. We searched 
also for the multidimensional maxima of the likelihood 



adopting a Downhill Simplex Method (Press et al. 
obtaining consistent results. 



1989) 



a 



0.5 




0.8 - 




0.0 0.2 0.4 0.6 0.8 1.0 



Fig. 2. — The Likelihood function of = Qm + f^A normal- 
ized to unity at the peak after marginalizing along the Qm — f^A 
direction. The dashed line is the cumulative likelihood. 



In Figure 2, we plot the likelihood of fl normalized to 
1 at the peak where, again, we have maximized along the 
i7_A/ — r^A direction. The likelihood shows a sharp peak 
near = 1 and this result is insensitive to the tradeoff 
between D,m and f2A- (see Figure 3 and explanation in 
following paragraphs). This is an extreme manifestation 
of the "cosmic degeneracy" problem (because we are fo- 
cusing on just the first peak): we are able to obtain robust 
constraints on Q without strong constraints on CIm and 
Q\ individually. 

Within the range of models we are considering, we find 
that 68% of integrated likelihood corresponds to 0.85 < 
n < 1.25 (0.65 < n < 1.45 at 95%). The best fit is a 
marginally closed model with VIcdm — 0.26, fis — 0.05, 
f^A = 0.75, ns = 0.95, h = 0.70, Cio = 0.9. An almost 
equivalent good fit is given by ^cdm = 0.39, Ob = 0.07, 
^A = 0.65, ns = 0.90, h = 0.55, do = 1.0. 

In Figure 3 (top panel) we estimate the likelihood of 
the data for a 20 x 20 grid in (flMj^A) by applying the 
maximization/marginalization algorithm described above. 
The effect of marginalizing is, as expected, to expand the 
contours along the ^constant lines but has little effect 
in the perpendicular direction and we are able to rule out 
a substantial region of parameter space. 

For ~ 1 models, the position of the peak is solely 
dependent on the angular-diameter distance, with a good 
approximation being ipeak oc fi^a; this approximation 
breaks down when ^Im — > where the early time inte- 
grated Sachs- Wolfe effec t becomes important and £peak is 
far more sensitive to ft (White & Scott 1996). This effect 
leads to a convergence of contour levels as r2j\/ ^ in 
Figure 3. 
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Fig. 3. — The likelihood contours in the flM,flA plane, eval- 
uated at the maxima of the remaining four "nuisance" param- 
etes. The top panel is from the Boomerang/NA data, the 
bottom panel is from Boomerang/NA-|-COBE. The contours 
correspond to 0.32, 0.05 and 0.01 of the peak value of the like- 
lihood. The small triangle indicates the best fit. The dashed 
line corresponds to the flat models. 



4. DISCUSSION 

In the previous section we have obtained a constraint 
on using only the Boomerang/NA data. These new 
results are consist ent with Lineweaver (1998), Tegmark 
(1999) and pK99| . However, the Boomerang/NA data 
on its own does not constrain the shape and amplitude 
of the power spectrum at ^ < 25 and limits our ability 
to independently determine the parameters ns, ils, h, fl^ 
and CiQ. We combine the Boomerang/NA data with the 
4-year COBE/DMR angular power spectrum to attempt 
to break this degeneracy. In Figure 3 (bottom panel) we 
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plot the likelihood contours, again maximized over the nui- 
sance parameters for the combined Boomerang/NA and 
COBE data. The inclusion of the COBE data does not 
greatly affect the constraints at high flj^j or the confidence 
levels on 17, but, as expected, it helps to close of the con- 
tours at low values of fl^j- The best fit model changes 
to have ncoM = 0.46, = 0.05, = 0.50, ns = 1.0, 
h = 0.70, Cio = 0.94. We find that for the likehhood to be 
greater than 0.32 of its peak v alue then > 0.2, again 
simil ar to the results of DK99|. 



of the Boomerang/NA test flight to constrain the curva- 
ture of the universe. Given that we have based our results 
on this data set alone, our results are completeley inde- 
pendent from previous analysis of the CMB. At the time 
of submission, this letter is also the first analysis of this 
kind to include closed models in the computation. 

We find strong evidence against an open universe: we 
find that 0.65 < fl < 1.45 at the 95% confidence level, sig- 
nificantly ruling out the current favourite open inflationary 
models for structure formation ( Lyth fc Stewart 1990|, Ra- 



from the l uminosity-distance meas u rements of high-z su - 



One can combine our constraints with those obtained tra & Peebles 1995, Bucher, Goldhaber & Turok 1995) 



pernovae (Perlmutter et al. 1998, Schmidt et al 199 



using the l-cr constraint from Perlmutter et al (1998 
f^M - 0.6f7A = -0.2 ± 0.1, we find 0.2 < < 0.45 
and 0.6 < f^A < 0.85. 

A few comments are in order about the robustness of our 
analysis. Firstly we have not truly marginalized over the 
nuisance parameters. However the constraints we obtain 
in this way are, if anything, more conservative. Secondly, 
although we are limiting ourselves to standard adiabatic 
models, a strong case can be made against the rival the- 
ory of topological defects: the presence of a fairly local- 
ized rise and fall in the data around £ of 200 indicates that 
the characteristic bro adening due to decoherence of the 
either cosmic strin gs (IContaldi, Hindmarsh fc Magueijc 
19990 or textures ( jPen, Seljak fc Turok 1997| ) is stron"gI^ 
disfavoured. 

Finally we have restricted ourselves to only four ex- 
tra nuisance parameters. Again we believe this does not 
affect our main result (our constraints on fi) although 
it may affect the low D,m constraints when we combine 
the Boomerang/N A data with COBE; the results from 
Tegmark (1998) and |DK99| lead us to beheve that the ef- 
fect will not greatly change our results. 

To summarize we have used the angular power spectrum 



Much tighter constraints will soon be placed on these and 
others cosmological parameters from future data sets, in- 
cluding data obtained during by the Antarctic LDB flight 
of Boomerang, which mapped over 1200 square degrees 
of the sky with 12' angular resolution and higher sensitiv- 
ity per pixel than the Boomerang/NA. 
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